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Abstract. 

Many astrophysical sources radiate via synchrotron emission from relativistic 
electrons. The electrons give off their kinetic energy as radiation and this radiative 
loss modifies the electron energy distribution. An analytical treatment of this problem 
is possible in asymptotic limits by employing the continuity equation. In this article, 
we are using a probabilistic approach to obtain the analytical results. The basic logic 
behind this approach is that any particle distribution can be viewed as a probability 
distribution after normalizing it (as is done frequently in statistical mechanics with 
ensembles containing very large number of particles). We are able to reproduce the 
established results from our novel approach. Same approach can be applied to other 
physics problems involving spatial or temporal evolution of distribution functions. 


1. Introduction 


Synchrotron radiation is emitted by relativistic charged particles circling magnetic held lines. 
Highly relativistic particles are abundant in our universe and the radiation detected from a range 
of astrophysical sources originate due to this process. Supernova remnants, pulsars. Gamma 
Ray Bursts and large-scale jets emanating from certain galaxies are some of the examples. In 
this article, we are describing the time evolution of a population of electrons emitting synchro¬ 
tron radiation. We used probability theory to arrive at the results conventionally derived using 
the continuity equation. These results are widely used in modeling astrophysical sources. Our 
approach is generic and can be used for similar physical problems. 

Total energy of a relativistic electron can be expressed as £ = ynieC^ where nieC^ is the 
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rest energy and y is the Lorentz factor ( -^1 - ) of an electron moving with velocity v. 

The synchrotron power depends on the energy E (equivalently lorentz factor y) and the 
larmor radius of the electron. As the electron radiates, its energy (equivalently y) reduces. Time 
evolution of the Lorentz factor of a synchrotron radiating electron of mass mg, charge e and 
initial lorentz factor yo in the presence of a magnetic field density B can be derived as follows 
(see appendix for the derivation) 


7(0 = 


yo 

1 + Ayot 


where cr, is the Thompson’s Cross section. 




( 1 ) 


Instead of monoenergetic electrons, realistic systems have electrons of a certain energy distri¬ 
bution. Most astrophysical systems emitting synchrotron radiation do not have sufficient particle 
density to achieve thermal equilibrium. Hence instead of a relativistic Maxwellian, the energy 
distribution is expected to be a ‘power-law’ of the form : 

n(y)dy = Key~’’dy where ym<7 <yu (2) 

where y^ and y^ are the minimum and maximum limits of the distribution respectively. Kg is an 
arbitrary constant depending on the physical parameters of a given system. This is also called the 
‘non-thermal’ distribution. This distribution does not remain the same over time due to synchro¬ 
tron loss suffered by each electron. However, the time evolution of the electron population not 
only depends on the energy loss but also on the injection of fresh particles into the system. Our 
aim is to obtain the distribution n'(y(0) at any given time t. 

We assume a constant magnetic field B and two types of injection : (i) One Shot injection 
and (ii) Continuous Injection. We assume that the fresh particles that are injected in the system 
always follow the non-thermal powerlaw given in equation 2. 


2. One Shot Injection 


One shot injection means that at t=0 we inject electrons in the system which follow the power-law 
distribution in equation 2. If we divide equation 2 by the total number of particles in the system 
N, it becomes a probability distribution ; 

p(7o)dyo = Ky^^’dyo where y„, < yo < oo (3) 

We have assumed that y^ oo. K = Kg/N and y^ is the minimum lorentz factor of the distribu¬ 
tion. p(yo) - is the probability that an electron has an initial Lorentz factor between yo and yo + 
dyo 


If we substitute yo as oo in equation 1, then y(f) will tend to Therefore, even if the 
maximum lorentz factor in infinite initially, it will to tend to some finite value say at a later 
time t. 
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The system at a later time t can be represented as, 

p(y)dy - Ky^^’dy where ym < y < yu (4) 

K - Kg IN', piy) - probability that a electron has lorentz factor between y and y + c/y at a given 
time t. 

Let us consider two random variables: 


1. yo which represents the Lorentz factor of an electron at time f = 0 

2. y which represents the Lorentz factor of an electron at some time t 


In order to fin d p(y), we can use the theorem for transformation of random variables from prob¬ 
ability theory ( MversllVe and Walnole 2007 ). 


Theorem 1. Suppose that X a continuous random variable with probability distribution f(x). 
Let y — u(x) define a one-one relationship correspondence between the values ofX and Y so that 
the equation y — u(x) can be uniquely solved for x to obtain x — w(y). Then, the probability 
distribution ofY is given is: 

giy) - f(w{y))\J\ (5) 

where J - dwidy and is called f/ie Jacobian of the transformation 


In our case x - yo, f{x) - p (yo), y - y and g{y) - piy), uix) is given by equation 2. So, 
using Theorem 1, we can write. 


piy) = p'iyoiy))^ ^k- ^ 


dy (l-Ayf)^"^^ 


or, 


piy)dy = K 


{\-Ayt)-P^^ 


dy 


where 


Ym <7 <yu 


(6) 


(7) 


At a later time t, minimum and maximum values of y have changed which can be found out 
by substituting initial minimum (y) and maximum (cxi) values in equation 2. 


Ym = 




1 + Ayt 




( 8 ) 


Note, that even if at f = 0, maximum possible Lorentz factor for an electron was infinite, after 
a time t, maximum possible Lorentz factor = y^. No electron can have a Lorentz factor greater 
than yc at a given time t. 
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lug (r) 


Figure 1. Initial and Time Evolved Particle Distribution for one shot injection given by equation 2 and 
equation 9 respectively 

Equation 7 gives us the probability distribution. For obtaining the particle distribution func¬ 
tion, we have to multiply the probability distribution function by N (total number of particles). 

n(y)dy^Ke- - ——^dy where Tm < 7 < Tc (9) 

(1 - Ayt) P 

where Kg - NK. So, the particle distribution function at time t for one shot injection is given 
by equation 9. Figure 1 shows both initial and time evolved distribution functions. The graph is 
plotted between y,„ and y^. 


3. Continuous Injection 


In the case of Continuous Injection, fresh electrons continuously get added and make up for the 
loss in energy space. Now, total number of particles will also vary. We can treat Continuous 
Injection as a series of one shot injections. For simplicity, let’s assume that at each instant of 
time we inject Nq number of particles and each injection follows equation 1 with same p,y,„ and 

r» = 


n(yo)dyo - K^y/dyo where y„, < yo < oo 


(10) 
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Ke is determined from the equation 



niyo)dyo - No 

(11) 

II 

(12) 


Now, again consider two random variables yo and y in exactly the same manner as described 
for the one shot injection case. Each electron in the system was injected during a particular injec¬ 
tion and each injection evolves independent of other injections. If we look at any one injection, 
it evolves in a manner similar to what we discussed in previous section for one shot injection. 
The overall evolution can be viewed as a summation or integration over all injections. If we are 
looking at time t we can assign two other variables to each electron f and f 2 , ty is the time when 
the electron was injected and f 2 is the time elapsed since its injection or the time for which the 
electron has been there in the system. Both t\ and f 2 varies from 0 to f such that h + t2- t. 

Define: 


piy, h) Probability that a random electron has a Lorentz factor y at time t and has been there in 
the system for time f 2 . 

p(y) Probability that a random electron has Lorentz factor y at time f. 

p(t 2 ) Probability that a random electron has been there in the system for time t 2 - 

piyih) Probability that a random electron has a Lorentz factor y given that it has been there in 
the system for time f 2 . 


Lrom definition of conditional probability distribution, we know that: 

piy, h) = p{t2)p{ylt2) 


As the injection is identical at every instant, p(t 2 ) is a uniform distribution: 


Pih) = 



if 0 < f2 f 

otherwise 


(13) 


(14) 


piyih) is the probability that a random electron has a Lorentz factor y given that it has been 
there in the system for time f 2 . If we know the initial distribution and time elapsed, then, the final 
distribution is given in a same way as for one shot injection case (Once we fix f 2 , we are talking 
about a particular injection, so the problem is identical to one shot injection case). Therefore, 

y^P 

piyih) - K - - where 0 < t2 < t (15) 

(1 -Ayt)~P^^ 
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where, K - Kg I No, putting equation 14 and 15 in equation 13; 




K y-P 

t (1 -AytyP^^ 


0<t2<t 


(16) 


So, far we have not imposed any condition on y. Now from definition of marginal distribution 
function we know that: 


piy) 




piy, t2)dt2 


(17) 


Substituting equation 16 in equation 17: 


piy) 


f 


K 


piy) 


K Y 


-p+\ 


A ip- \)t 


1-(1 


rd 

+ 

1 

(18) 

1 

^■72 

1 

(19) 


yc and are defined by equation 8. 

Depending upon yc and y^, we can have two cases: slow cooling (y^ < y^) and fast cooling 
iym > yc)- 


The loss becomes considerable only for those electrons for which the radiative time scale is 
less than the age of the system ityad < 0 because only those electrons will dissipate energy at a 
faster rate than the age of the system. Radiative time scale is given as: 


E 

" dEjdt 


( 20 ) 


So, loss becomes considerable for those electrons for which 

E 

dEjdt ~ 

From equation 1, we can easily derive; 

dt 1 + Ayof 


( 21 ) 


( 22 ) 


Using above equation and the relation E - ynieC^, we can easily derive the following relation: 


dE 

dt 


diyiHcC^) 


dy 2 

— nieC 

dt 


■■ Ay^nicC^ 


dt 


(23) 
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Substituting above equation in equation 21 we get, 

ynieC^ 

Ay^tfieC^ 

Hence, Equation 21 is equivalent to equation 24. Both are just different representations of 
the same concept. Both implies that the energy loss is greater for electrons having a lorentz factor 
greater than jc, we expect the power law to be steeper in the range y > yc- 

Therefore, if y^ < yc, only a small fraction of electrons in the system will be cooling dom¬ 
inated. Hence we refer it to as slow cooling. However, if y^ > yc, a large fraction of electrons 
is cooling dominated and hence we refer it to as fast cooling. The electron distribution at a given 
time is different for both cases. Hence, we have discussed them separately in coming subsections. 


r > 


At 


(24) 


3.1 Slow Cooling 

The expression in equation 19 is valid only for ym < 7 < 7c- From equation 8,we know that the 
maximum attainable lorentz factor for a electron at a time t is given by yc = 1 lAt. A continuous 
injection can be seen as a superimposition of various sequentially injected one shot injections. 
Hence, we can apply the results derived for one shot injection on a particular sub-injection of the 
continuous injection. 

Consider a continuous injection started at f = 0 and current time as t. Now, take a one shot 
injection injected at time ti into the system. Then, any electron which was a part of this injection 
has been inside the system for t 2 - t - t\ seconds. Therefore, the maximum attainable lorentz 
factor for such a electron is 1 /Af 2 . 

Suppose we want to calculate p{y) for some y > y^ = 1 /At. Only those electrons can attain 
this value whose maximum attainable lorentz factor is greater than the value of y at which we are 
calculating p{y) i.e. electrons for which 


— > y 

At2 

(25) 

1 

t2 < — 

Ay 

(26) 


Therefore, for calculating p{y) for y > y^, upper limit of the integral in the equation 18 will 
change to 1 /Ay instead of t. 


Hence, for y > yc 


piy) 


_ K 

~ Jo ”(1^ 


(1 -Aytyp*^ 


dt2 


(27) 
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piy) = 


K 

A(p - l)t 




7>rc 


(28) 


Equation 19 and 28 combined gives p(y) for < y < oo, but there will be some electrons 
which will have Lorentz factor below also (electrons having initial Lorentz factor just above 
ym will lose energy to go below y^. For these electrons also, we will have a different expression. 
The minimum Lorentz factor attainable at time t is y^, which is given by given by equation 8. 
But this is the minimum for the very first injection injected at f = 0. For electrons belonging to 
any later injection (which were injected at time fi > 0 and has been there for time t 2 < f), the 
minimum attainable Lorentz factor is ym/(l +Aymh)- So, if we are looking at a y < y^, only those 
electrons have a non-zero probability of having this Lorentz factor which were injected such that 
the minimum attainable Lorentz factor at present time is less than y. 


Jn. 


1 -H Ajmtl 

y-y 


<y 


h > 


yyn. 


(29) 

(30) 


So, for calculating p(j) for y < y^, lower limit of the integral in the equation 18 will change 
to (y - ym)lyym instead of t. So, for y < y;„ 


p(y) = 


r-/ 

yym 


K y-'’ 

- 

t (l-Ayt)-P^^ 


p(y) = 


KyAP+l) 

A(p - l)t 




y <ym 


(31) 

(32) 


3.2 Fast Cooling 


The main logic remains the same for fast cooling also, equation 16 and 17 are always valid. If you 
are looking in the region where y > y^, then, upper limit in the integral in equation 18 should be 
as given by equation 24 instead of t. Similarly if you are looking in the region where y < y^, the 
lower limit of the integral in equation 18 should not be zero. But, as given by equation 30. Now, 
in case of fast cooling, if you are looking at the region yc < y < ym, only those electrons have 
a non-zero probability of having this Lorentz factor which were injected such that the minimum 
attainable Lorentz factor at present time is less than y and the maximum attainable Lorentz factor 
at present time is greater than y i.e. injections which satisfies both equation 26 and 30. So, both 
the limits of the integral in equation 18 will change according equations 24 and 30. 


So, for fast cooling and ym <y < yc. 


p(y) 


X 


^ K 


t (1 - Ayf) 


dt2 


(33) 
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piy) = 


K y-2 


7m <7 <yc 


(34) 


For y > y„, only upper limit changes because y > 7c but not < so the expression for p(y) 
is exactly same as equation 26. 


Similarly, for y < yc, only the lower limit changes because y < ym but not y > yc, so the 
expression is exactly same as equation 30. A summary of formulas for p(y) in different regimes 
for both slow and fast cooling case is given below: 


Slow Cooling 


p(y) = 

p(y) = 


Ky 


<P+i) 


A(p - l)t 


p-i 


K y-P-^^ 

A(p-l)t 


1-(1 


p(y) = 


K 


A(p - l)f 



yc 


.y-(P+l) 


P-n 


ym 


y>yc 


y <ym 


<y <yc 


(35) 


(36) 


(37) 


Fast Cooling 


piy) = 


Ky 


.-(p+i) 


A{p - \)t 


p-i 


-11-^ 


p-i 


y<yrr, 


(38) 


piy) 


piy) = 


K y-2 

I(^^ 

K 


Aip - \)t 


y 


(/>+!) 


ym<y <yc 


y>yc 


(39) 


(40) 


These are probability distributions, for obtaining particle distribution functions, multiply all 
of these equations by A = NqI (total number of particles). It will simply change Kjt in all the 
above equations to Kg. 

Not X — = NqK - A),(mentioned earlier) (41) 
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The particle distribution functions are given below; 
Slow Cooling 


n(y) 


A{p - 1) 


p-i 


p-i 


7 <yir. 


n(y) 


Ke 

A {p-D 


y p-i 

1 -(!--) 

7c 


ym <7 ^ 7i 




7 > 7c 


(42) 


(43) 


(44) 


Fast Cooling 


n(y) 


KeyAP^^^ 

A(p - 1) 


xf-b-i''’" 


7 < 7n : 


(45) 


n(r) = 


Ke y-" 

Mp-DyP^^^ 


7m <7 < 7c 


(46) 


K 


£_^,-(P+l) 


”M=S(7TTjX'"' r>r. 


(47) 


Figure 2 and Figure 3 shows the variation of particle distribution as a function of y for the 
slow cooling case and fast cooling case respectively in logarithmic scale. 


4. Standard Result 


The standard result ( Wallace & 1979lChapter 3h obtained by solving the continuity equation says 
that electron distributions are multi-powerlaws with power indexes defined as follows: 


• p for ym<y < 7c and p+1 for yc<7 

• -2 for ym < 7 < 7c and p+1 for yc < 7 
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Slo.s Ceding 



Figure 2. Particle Distribution Function for Slow Cooling case given by equations 42-44 

Analytical Solutions exists only in approximate regimes and full result requires numerical calcu¬ 
lations. 


5. Conclusion 


We can see that the results obtained by this theoretical probabilistic model are similar to the ones 
obtained by solving the continuity equation, p + I power-law index part is exactly same for both 
slow cooling and fast cooling. Also, 2 power-law index part is same for fast cooling and for 
much smaller values of y, even equation 43 approximates to a power-law index p. However, 
probabilistic model opens a new dimension, as it is intermediate to the analytical solutions like 
Tucker 1978. The probabilistic approach to this problem has not been attempted so far to the 
best of author’s knowledge. Such an approach is possible in similar problems which involves the 
evolution of a distribution. 

One striking feature observed in these expressions is the steep rising part (y < y„, in Fig. 
2 and 3) which is normally not discussed much. The reason for this is that y„, and are very 
close and the fraction of electrons is negligible. But, the model discussed here clearly suggests 
that there will be some electrons between y,„ and y,,, which shows a rising particle distribution as 
we go towards higher y. However, this behavior is confined to a very narrow range of y. If the 
magnetic field is very high then this sharp rise becomes significant as the difference between y,„ 
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Fast Cooling 



Figure 3. Particle Distribution Function for Fast Cooling Case given by equations 45-47 


and will be more and a larger fraction of electrons will lie in this region. Theoretically this 
behavior should be present in electron distributions evolving under a magnetic field. 

These results were obtained assuming a constant magnetic field for simplicity. But, the gen¬ 
eral concept is very simple and logical. Even for a time varying magnetic fields we can apply 
exactly the same procedure. The final expressions will change depending on the functional form 
of B but, all the integrals/integral equations will be same. For e.g. equations 13, 14 and 17 are not 
going to change under any circumstances. Also, the logic behind changing the limits in different 
regimes of y for both slow cooling and fast cooling (equation 24 and 28) is also the same however, 
the final conditions change depending upon what are the minimum and maximum values of y at 
a later time t. 

The main logic behind the approach is that any particle distribution problem can be modeled 
as a probability problem because both are similar mathematically (just a difference of normal¬ 
ization). In a generic sense, if we have a particle distribution and it is varying with respect to 
any variable (spatial or temporal), we can model it as a probability problem by normalizing the 
particle distribution. 
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A. Appendix material 


Synchrot ron Power emitted by an ultra-rel ativistic electron under a constant magnetic field is 
given by( Rvbicki & Lightman IQTQjfDg 169 ) 


where t/a is magnetic energy density = fracB^2>n. Also, we know that: 

dE d , 9, 

P = = -—{ynieC ) 

dt dt 

Substituting Equation 49 in Equation 48: 


idy 4 ^ 

= = —fi. 


' 2 ,, 2 ^ _ ^tC 2,2 


Separating the variable y and f, and then integrating we get: 

cr,B^ 


r. 

Jyo 7^ Jo 


67mieC 


dt 


1 1 


t - -At 


Jo j 6mneC 

r - To = -yjoAt 

ro 


1 + Ajot 


A = 


6nmpC 


(48) 

(49) 

(50) 

(51) 

(52) 

(53) 


r (0 = 


where 


(54) 
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